Mon. Not. R. Astron. Soc. 000, [[H5] (2010) 



Printed 2 September 2010 



(MN m$K style file v2.2) 



A model for nulling and mode changing in pulsars 



A. N. Timokhin 1 ' 2 * 



O 

Oh' 

c/3 



X 

6 

o3 



> 
in 

On 

on 
o 



1 Astronomy Department, University of California at Berkeley, 601 Campbell Hall, Berkeley, CA 94720, USA 
2 Sternberg Astronomical Institute, Universitetskij pr. 13, Moscow 119992, Russia 



Received ; in original form 



X 



ABSTRACT 

We propose that in some pulsars the magnetosphere has different states with different 
geometries or/and different distributions of currents, it occasionally switches between 
them. These states have different spindown rates and emission beams, in some of the 
states no radioemission is produced at all. Switching into a different state manifests 
as a mode change when we see different parts of the emission beam or the beams in 
different states have significantly different geometries, it manifests as nulling when we 
cither miss the new beam or no radioemission is generated in the new state. We show 
that modest variations in the beam shape can be accompanied by large variations 
in the pulsar spindown rate W - the dependence of W on the opening angle of the 
emission beam a can be as strong as W oc a 4 . We speculate about physical mechanisms 
which may cause reconfiguration of the magnetosphere. 
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1 INTRODUCTION 

Rotationally powered pulsars are observed mostly in radio 
band because of the very high sensitivity of radio telescopes. 
Pulsar radioemission shows very rich set of phenomena. Ex- 
amples of such phenomena are mode changing and nulling, 
which are observed in many pulsars. There are hints that 
mode changing and n ulling may be ma nifestations of the 
same phenomenon (e.g. IWang et al.ll2007l ). It is not clear yet 
whether nulling and/or mode changes are 'uniform' phenom- 
ena, i.e. whether there is a single class of reasons why pulsars 
null or/and change modes. There are nulling pulsars which 
are quiet for only a few periods and there are pulsars which 
spend very long time in each of the state. There are pulsar 
for which the emission intensity during nulls drops at least 
by several orders of magnitude or switches off completely so 
that they cannot be detected even in a deep search, while 
for some pulsars nulls seems to be just mode(s) wi th low 
emission intensity (e.g. IWang et alj|2007l : lBiggjll992h . 

It was not clear for a long time whether nulling is due to 
'micro-physics' of the radioemission mechanism or changes 
in the pulsar magnetosphere as a whole play a role. The 
discovery of two nulling pulsars PSR B1931+24 and PSR 
J1832+ 0029 with different spindown ra tes in ON and OFF 
phases IjKramer et alj|200q ; lLvndl2009l ) provides strong ev- 
idence that at least some subclass of nulling may be re- 
lated to the global properties of pulsars and not merely a 
property of radioemission mechanism(s) alone. In this short 
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paper we will focus on pulsars where global processes in- 
volving changes in the whole magnetosphere seems to be at 
play. We propose a qualitative model which could explain 
nulling and mode changing as manifestation of the same 
phenomenon and account for the observed variation of the 
spindown rate. 



2 THE MODEL 

The radioemission of a pulsar is negligible in the overall en- 
ergy budget, usually constituting a very tiny fraction of the 
pulsar spindown rate, less then 10 -3 . Most of the energy 
flux is carried away by the relativistic pulsar wind and does 
not reveal itself in the emission. Therefore, merely the fact 
whether the radioemission is produced or for some reason is 
switched off or strongly suppressed cannot change the pul- 
sar spindown such that it would be measurable. However, 
there are at least two pulsars - PSR B1931+24 and PSR 
J1832+0029 - which exhibit nulling and at the same time 
their spind own rates are substantially less when they are 
non-visible (|Kramer et al ] |200rj ; |Lvnefl2009l ). In the case of 
PSR B1931+24 the difference in the spindown rate between 
the ON and OFF modes is about 50%. Nulling in these pul- 
sars takes an extreme form, they spend a very long time 
in both state (days), whereas many other nulling p ulsars 
switch off for much sho rter times - hours or less (e.g. iBiggsl 
1 19921 : IWang et ai1l2007l ). Taking into account the smallness 
of the pulsar period derivatives, P ~ 10~ 15 , it is very hard 
to detect variations of spindown rate between modes lasting 
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only hundreds or thousands of pulsar periods, even if these 
variations are large. We may speculate that other nulling 
pulsars, or at least some of them, can have non-negligible 
variations of the spindown rate too, but they are difficult to 
measure and are not detected. 

The behavior of PSR B1931+24 and PSR J1832+0029 
hints that nulling, at least in some pulsars, is a manifes- 
tation of some global process which changes the spindown 
rate of pulsar. This global process should involve changes in 
the structure of the pulsar magnetosph ere. The most pop- 
ular r adiopulsar model, introduced by iGoldreich Sz Julian! 
(1 19691 ). implies that the pulsar magnetosphere is filled with 
plasma, and accelerating electric field is present only in ge- 
ometrically very small regions - in the pulsar polar cap 
and in regions around and inside the current sheet carry- 
ing the bulk of the return current. The rest of the magneto- 
sphere should be force-free. In the force-free magnetosphere 
the energy losses of pulsar are determined by the configu- 
ration of the magnetosphere. Within this configuration we 
understand both the geometrical relation between the open 
and closed magnetic field line zones as well as the distri- 
bution of the current density in the open field line zone. 
In other words, in the force-free model changes in the pul- 
sar spindown rate are possible only if the configuration of 
the magnetosphere changes. Recently the force- free pulsar 
magnetosphere model has been studied in great detail (e.g . 
Contopoulos et al.ll 19991; lGruzinov | l2005| ; IContopoulosll2005l ; 
Timokhinll2006l . l2007al ; ISpitkovskvboOrj) . For an aligned ro- 
tator it was shown that there exists: i) stationary magneto- 
spheric configuratio ns with different sizes of the closed mag- 
netic field line zone (|Timokhinll2006l ); and ii) configurations 
with diffe rent current density distributions in the open field 
line zone (|Contopoulosll2005l ; lTimokhinll2007al lbh . These con- 
figurations have different spindown rates. 

Having the above described facts, hints and specula- 
tions in mind we make the following conjecture. Magneto- 
spheres of some pulsars at some stages of their evolution 
can have sets of quasi-stable states with different sizes of 
the closed field line zone or different current density distri- 
butions in the open field line zone, or both. Occasionally a 
magnetosphere switches between these states. 

Currently there is no reliable theory which can explain 
pulsar radioemission. The radioemission is coherent and it 
seems that conditions for generation of coherent emission in 
pulsar magnetosphere may be not always fulfilled - at least 
there are pulsars with periods and magnetic fields typical 
for radioemitting pulsars which are v isible in gamma-ray s 
but show no detectable radioemission l|Abdo fc et al.| [2010). 
If, at least for some range of pulsar parameters, the condi- 
tions for generation of coherent emission are very sensitive to 
the properties of the magnetospheric plasma - like current 
and plasma densities or/and their gradients - then it can 
be that in some of the magnetospheric states the coherent 
emission is not generated at all or is generated at different 
places. Each state then has unique radioemission properties 
and switching between different magnetospheric states will 
manifest as mode changing or nulls. 

Another reason for changes of the pulsar mean profile 
can be due to self-similar changes of the emission beam when 
the radioemission zone shrinks or expands, i.e. radioemission 
is generated in the same way in all states but the beam 
scales self-similarly with the size of the emission zone. Then 



the emission beam in some states is wider than in another. 
Depending on the orientation of our line of sight relative 
to the emission beam we can see either mode change or 
null depending on whether our line of sight still crosses the 
narrower emission beam or not. 

We note that in the proposed model the magnetosphere 
is always filled with plasma and qualitatively works in the 
same way, i.e. particles are accelerated, gamma-rays are 
emitted etc. But changes in the beam pattern due to current 
redistribution or/and shrinking of the corotating zone are al- 
ways accompanied by changes in the spindown rate. Below 
we try to quantify the proposed model. We defer specula- 
tions about reasons causing such behavior of the magneto- 
sphere to Discussion. 

We do not know what the conditions for coherent emis- 
sion of radio waves in the pulsar magnetosphere are. Hence, 
we can not make any quantitative statements about what 
changes in the magnetospheric configuration can cause fail- 
ure of the emission mechanism or what is the necessary 
condition for large changes in the beam pattern when the 
magnetospheric configuration changes. However, we can esti- 
mate how the emission beam changes when the emission pat- 
terns in all configurations are similar - when nulling/mode 
changes are due to pure geometrical effects. Even if con- 
traction/expansion of the emission beam takes place only in 
some nulling/ mode changing pulsars, detection of a correla- 
tion between the width of the mean profile and the spindown 
rate in a mode-changing pulsar could be a strong argument 
in favor of the proposed model. On the other hand, for pul- 
sars where nulling/mode changes cannot be (fully) explained 
by contraction/expansion of the emission beam, such esti- 
mates could provide useful limits on variation of the spin- 
down rate. 

For simplicity we consider only the case of an aligned 
rotator with a dipole magnetic field. The main point we 
wish to demonstrate with these crude estimations is that 
large variations of the spindown rate can be accompanied 
by modest changes in the geometry of the radio beam, or, 
reversing the statement, if changes of the beam width due 
to reconfiguration of the force-free magnetosphere are ob- 
served, then changes in the spindown rate would be large 
and could be detected. 



2.1 Shrinking-expanding corotating zone 

Let us assume that a pulsar magnetosphere has different 
states with different sizes of the corotating zone, the zone 
with closed magnetic field lines. Let us consider two states: 
state A with the corotating zone smaller than that in the 
state B (see Fig. [TJ. The size of the polar cap, limited by 
the last closed field line is smaller in state B. The current 
density distribution in the polar cap of p ulsar for such stat es 
do not differ dramatically (see Fig. (5) in lTimokh in (2006)). 
It is naturally to assume that all processes in regions not far 
from the neutron star (NS) for small changes in the size of 
the polar cap behave similarly, i.e. the beaming pattern in 
different configurations can be approximately described by 
a shrinking/contraction proportional to the changes in the 
size of the polar cap. 

If radio emission is generated close to the NS, then the 
emission beam in state A has a larger opening angle. When 
the pulsar switches from state A into state B, depending on 



Nulling and mode changing in pulsars 3 



Mode Change 




Figure 1. Schematic view of pulsar magnetosphere when the 
closed field line zone changes its size. Radioemission beams are 
shown by the blue lines, the last closed field lines - by the black 
lines. The lines of sight for which pulsar will demonstrate mode 
change and null are shown by magenta and red arrows respec- 
tively. 



our line of sight, we see either different parts of the emission 
cone, what might result in a mode change, or we miss the 
emission cone entirely and conclude that pulsar is in the 
'null' state. Configurations with a smaller corotating zone 
have more open magnetic field lines and, hence, have higher 
energy flux to infinity. Therefore, during the ON state the 
spindown rate of pulsar is higher. 

If emission is directed along magnetic field lines, then 
the opening angle of the emission cone a is 



a = 9 + arctan (B e /B r ) ~ 3 9/2 oc 8 P 



(1) 



where 9 is the colatitude of the emission zone, Be , B r are 
components of the magnetic field in spherical coordinates, 
and 6 pc is the colatitude of the polar cap boundary. The last 
step in eq. JT]) comes from our assumption about similarities 
of emission zones in different configurations. Energy losses 
of the aligned rotat or depend on th e size of the corotating 
zone as (eq. (62) in lTimokhinl lj200rj ^ 



W~ W md 



(— V 



(2) 



where Ry is the actual size of the corotating zone, Rlc is 
the size of the light cylinder. W m d are the magnetodipolar 
energy losses W m d = -Bq-Rns^ns/4c 3 , where B is the mag- 
netic field in the polar cap, i? NS is the NS radius, f2 N s is the 
pulsar angular velocit y, c is the speed of light. For dipolar 
magnetic field 9 pc ~ TyR NS /Ry. Expressing Ry through 9 pc 
we get 



W ~ W u 



(RlcY 



9 4 



V Rms J 

Taking into account relation $Q we finally get 
W oca 4 . 



(3) 



(4) 



We see that reconfiguration of the magnetosphere can indeed 
cause much smaller changes in the beam opening angle than 
in the spindown rate, and so small changes in the emission 
geometry can be accompanied by measurable changes of the 
spindown rate. 



2.2 Changing current density distribution 

Let us now consider the case, when the polar cap size is 
nearly the same (Ry does not change), but instead there are 
magnetospheric states with different current density distri- 
butions in the open field li ne zone . A m odel of such type was 
qualitatively discussed by Arons (1983) in relation to mode 
changing in PSR B0809+74. 

The force-free magnetosphere starts at the top of the 
cascade zone in the polar cap, above the pair formation 
front, where there is enough plasma to short out the ac- 
celerating electric field. The angular velocity in the magne- 
tosphere is constant along magnetic field lines. Each open 
magnetic field line rotates with the angular velocity Qf it 
has at the top of the cascade zone in the polar cap. Q,f de- 
pends on the distribution of the potential drop in the cascade 
zone Vj| - the potential drop between the base of the force- 
free magnetosphere and the NS surface. In the force-free 
magnetosphere the current density distribution is set by the 
requirement that the open magnetic field lines pass smoothly 
through the light cylinder. For magnetospheric states with 
different distributions of the accelerating potential Vj| the 
distributions of Q,f will be also different, and so will be the 
positions and the shapes of the light cylinder. Hence, such 
states will have different current density distributions in the 
open field line zone as well. 

Pulsa r energ y losses are given by (cf. eq. (40) in 
iTimokhinl (|2007al )) 



W ■ 



mi 

2nc 



(5) 



where J is the total current flowing through the polar cap, 



Bo7rr pc is the total magnetic flux in the open field 



line zone, r pc is the polar cap radius. The dependence of 
J on Q.f could be approximated by that dependence in 
the magnetosphere of a r otating split-monopole J oc Q.f 
|Blandford fc Znaield[l9"77l y And so we have for the energy 
losses 



W oc r 



(6) 



the spindown rate is proportional to the square of the total 
electric current flowing through the polar cap, a dependence 
typical for power dissipated in an electric circuit. 

To change the spindown rate the total current must 
change; small local changes in the current density which do 
not change the total current cannot account for changes in 
W . As a toy model let us consider two states A and B, where 
the shape of the current distribution remains the same but 
the current density j in the state B is less than that in the 
state A by some factor r\ slightly less than 1, 



Jb = VJA ■ 



(7) 



The total current changes then as Jb = r] J a , and for changes 
in the energy losses we have 



W oc t) 



(8) 



There is no reliable model for pulsar radioemission. We 
can only speculate what sets the pattern for the emission 
beam. Let us assume that in state B the emission zone shifts 
to the field line where the current density is equal to the 
current density which flowed through the emission zone in 
state A, i.e. 
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J B (fe) = JA(0A) 



(9) 



This toy model results in changes of the beam geometry 
similar to that shown in Fig. [1] 

The current density in the polar cap of the aligned ro- 
tator with negligible differen tial rotation o f the open mag- 
netic field lines is close to the iMichell (|l973T ) current density 
|Timokhin| [2006) 



3 =3G3 (1 



(10) 



where j G j is the Goldreich- Julian current density. Using re- 
lations and ([9]) with the current density given by eq. (|10[) . 
we get the following expression for variation of the emission 
zone colatitude in the case of small changes in the current 
density 



60 



7«& 



(11) 
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(12) 



From eq. ((8| we have SW/W = 2(77 — 1), and small variation 
of the energy losses are related to variation of the emission 
zone colatitude as 

SW 
W 

where ( = {6 2 / 9 pc ) / (1 ~ 6 2 / 6p C ) . If the emission zone is close 
to the edge of the polar cap, 8 > 9 pc / v2 ~ 0.71 8 pc , then £ > 
1. Taking into account the relation between the colatitude 
of the emission zone and the beam opening angle fU we get 
the power law 

W oc q 4C , (13) 

which is valid for small variations of W and a. We see that in 
this case changes in the spindown rate could be even larger 
than in the case of a changing corotating zone discussed 
before. Although the model considered here is very simple, it 
demonstrates nevertheless that the dependence of the energy 
losses on the beaming angle can be rather strong. 



3 DISCUSSION 

Let us summarize the reasoning behind the proposed model. 
It seems that the force-free model agrees well with observa- 
tions: pulses are narrow, pointing to smallness of the regions 
with accelerating electric field; as it is inferred from obser- 
vations of pulsar wind nebulae, there should be more than 
enough plasma in the magnetosphere to screen the accel- 
erating electric field. If pulsar magnetosphere is essentially 
force- free, then changes of the spindown rate are possible 
only if the configuration of the magnetosphere changes. So, 
in the case of PSR B1931+24 and PSR J1832+0029 we must 
conclude that we have to do with changes in the magneto- 
sphere configuration, or the force-free model fails. Different 
configurations should have different emission-beam geome- 
tries. 

We do not suggest that all nulling pulsars switch be- 
tween states with different emission beams, neither that 
changes in the emission beam could be always attributed 
to the pure geometrical effect of shrinking/expansion; the 
latter can work only for pulsars where the line of sight 
crosses the outer parts of the emission beam. In the case 
of PSR B1931+24 the line of sight seems to make central 
traverse of the beam (Rankin, private communication), and 



so the simple geometrical models of Sec. l2.lll2.2l do not work 
for it. In the frame of the proposed model nulls in this pul- 
sar should be then attributed to cessation of the radioemis- 
sion due to changed parameters of magnetospheric plasma 
in the OFF state. But if pulsar magnetosphere can switch 
between different states, shrinking/expansion of the emis- 
sion beam can take place in some pulsars - at least in the 
mode changing pulsar PSR B0943+10 with tangential sight- 
line tra jectory the modes have differe nt widths of the mean 
profile l|Rankin fc SulevmanovallioOtil ). The point we wish to 
emphasize here is that reconfiguration of the magnetosphere 
can result in modest changes in the emission beam and in 
substantial changes in the spindown rate at the same time, 
i.e. the emission beams in different states can look similar 
while the energy loss rate changes a lot. 

An important consequence of the proposed model is 
that there must be mode changing pulsars showing substan- 
tial variations of spindown rate, so to say, the mode-changing 
"twins" of PSR B1931+24 and PSR J1832+0029. If the 
magnetospheres can have different states, switching between 
the states should manifest as mode changing or nulls. The 
spindown rate variation in the above mentioned nulling pul- 
sars strongly suggests that such states exist. Hence, there 
also must be pulsars where states with different spindown 
rates manifest as different emission modes. The proposed 
model also implies that nulling and mode changing pulsars 
should be among the pulsars with high timing noise. 

Let us now discuss what could be the underlying phys- 
ical mechanisms for the proposed model. It is not a ques- 
tion whether for a given pulsar period and magnetic field 
strength different configurations of the magnetosphere are 
admitted or not. At least for the aligned rotator it has been 
shown explicitly that such configurations exist, and we see 
no reasons why this can not be the case for the inclined ro- 
tator. The real question is why the magnetosphere can have 
a set of meta-stable configurations, i.e. configurations where 
the magnetosphere can stay for times much longer than the 
pulsar period. 

Although the properties of the force- free magnetosphere 
are quite well known, we cannot say what is the config- 
uration of a pulsar magnetosphere. The reason for this is 
that the configuration of the magnetosphere depends on the 
"boundary conditions" - physics of the polar cap cascade 
zone and physics of magnetic reconnection in the current 
sheet, especially at the so-called "Y"-point where the current 
sheet of the outer magnetosphere merges with that along the 
boundary between open and closed magnetic field lines. The 
distribution of the potential drop in the polar cap cascade 
zone sets the angular velocity of the open magnet ic field 
lines fi F (e.g. IContopoulosl 120051 ; iTimokhinl kooVal lbl). the 
reconnection rate at "Y" -point might influence the size of 
the of the corotating zone (e .g. IContopoulos fc Spitkovskvl 
120061 ; iBucciantini et alj|200a ). The physics of both of these 
regions is poorly understood, so we can only speculate about 
what properties the whole system consisting of the force-free 
magnetosphere + the cascade zone + the current sheet can 
have. Here we suggest one possible scenario for existence of 
several meta-stable magnetosphere configurations. 

The total energy of the magnetosphere of the aligned 
rotator with Of = ^ns monotonically depends on the 
size of t he corotat i ng zo ne, it decreases with the increase 
of R Y l|Timokhinl l2006h . The total energy of the split- 
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monopole magnetosphere with different current densities 
mon otonically increase s with the increase of the total cur- 
rent (|Timo khin 2007a). It can be that the combining effect 
of different current density distributions and different sizes 
of the corotating zone, especially in the inclined rotatoiQ, re- 
sults in a set of meta-stable magnetospheric configurations, 
in the sense that such configurations represent local min- 
ima of the total energy - any small deviations of the current 
density and/or the size of the corotating zone from that in a 
given "minimal" configuration would result in a larger value 
of the total energy. 

Changes in the current density distribution require vari- 
ations of the potential drop in the polar cap. These varia- 
tions should be of the order of several percent of the poten- 
tial drop across the polar cap in order to produce noticeable 
changes in Qf (Timokhin 2007a). In young pulsars, where a 
very small fraction of the vacuum potential drop would be 
enough to short out the accelerating electric field, these vari- 
ations should be small. With aging of a pulsar, the freedom 
in the current density distribution will increase, thus making 
existence of multiple metastable configurations of the mag- 
netosphere possible. In the frame of the proposed model this 
might be the reason why mode changes and nulling are seen 
in older pulsars. 

The pulsar magnetosphere represents a case of a highly 
non-linear system with a complicated means of adjustment 
between its subcomponents. For example, the polar cap cas- 
cade zone should adjust the current density flowing through 
it to that required by the magnetosphere. It is not clear 
now how this adjustment works. The characteristic timescale 
of the cascade zone is of the order of microseconds, while 
the magnetospheric characteristic timescale is of the order 
of the pulsar rotation period. Pulsar magnetosphere being 
highly non-linear system with very different characteristic 
time scales might exhibit rich dynamical behavior. In par- 
ticular, some of these meta-stable configurations may be a 
sort of strange attractors, where a system spends substan- 
tial time and then suddenly changes to a different state. 
The time intervals between these changes may be very large, 
much larger than any characteristic time scale of the system 
and transitions can happen quasi-periodically. This might 
account for quasi-periodicity in changes between ON and 
OFF states in PSR B1931+24. 

Rotating radio transients (RRATs) - recently discov- 
ered transient radi o sources - are though t to be rotating 
neutron stars too (|McLaughlin et al] 12006 V The discovery 
of pulsar PSR J0941-39 which switches between a RRAT- 
like isolated burst-emitting mode and a mode with emission 
typical for a nulling pulsar might imply that RRATs are 
an ex treme case of nulling pulsars |Burke-Spolaor fe Bailed 
|2010T ). In the frame of the proposed model it means that 
RRATs are nulling pulsars which spend most of their time 
in a state(s) where our line of sight does not cross the beam 

1 In the time-dependent numerical simulation s of th e in- 
clined rotator mag netosphere by ISpitkovskvl feooeft and 
Kalapotharakos & Contopoulos (2009) constant flp = £7 N s and 
fast reconnection rate at Y-point were implicitly assumed. There- 
fore, the fact that the only stable configuration which was found 
in these simulations was the configuration with Ry = RhC does 
not exclude possibility that other stable configurations with dif- 
ferent Q,p and Ry exist. 



or no radioemission is generated. The existence of such pul- 
sars seems to be a natural consequence of our model. 

As a direct hint that the pulsar magnetosphere can in- 
deed evolve on timescales much larger than the rotational 
period, we can consi der the evolution of subpulse drift rate 
in PSR B0943+10 l|Rankin fe Sulevmanovaf l2006h . In this 
mode-changing pulsar the subpulse drift rate evolves after 
the offset of the "B" mode with a characteristic timescale of 
some ~ 4000 rotational periods. This time is larger that 
any characteristic timescale the force-free magnetosphere 
can have. It might be that in this case we directly observe 
manifestations of some non-linear process(es) involving ad- 
justment between the magnetosphere, the polar cap cascades 
and/or the current sheet. If the pulsar magnetosphere is in- 
deed a highly non-linear system, then the existence of differ- 
ent metastable states seems to be a quite natural assump- 
tion. 
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